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Abstract: Large- N matrix models coupled via multitrace operators are used to define, 
via appropriate double-scaling limits, solvable models of interacting multi-string theories. 
It is shown that although such theories are non-local at the world-sheet level they have a 
simple description of the spacetime physics. Such theories share the main characteristics 
of similarly coupled higher-dimensional CFTs. An interpretation has been given in the 
past of similar continuum limits in terms of Liouville interactions that violate the Seiberg 
bound. We provide a novel interpretation of this relation which agrees with the current 
understanding of Liouville theory and analogous observations in the AdS/CFT correspon- 
dence. 
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1. Introduction 



1.1 Setup and questions 

In the context of the AdS/CFT correspondence, it has been argued fH that the holo- 
graphic dual of the product of k conformal field theories (CFTs) in d dimensions, deformed 
by a set of multi-trace interactions that couple the CFTs together, is a theory of quantum 
(multi)-gravity (or better a multi-string theory) on a union of k AdSd+i spaces with a 
formally common boundary - the boundary being isomorphic to the space where the dual 
product CFT lives. Assuming that each CFT has an independent gauge group G{ and that 
there are no fields charged under more than one group, multi-trace interactions are the 
only way to couple the CFTs. We will assume in this paper that we are working in the 
large- iV limit. 

From the quantum gravity point of view this is an interesting setup for the following 
reasons: 

(i) This is a non-trivial example of an interacting multi-graviton theory with a UV 
completion (the completion provided by the dual CFT). 1/N effects generate a non- 
vanishing potential for the gravitons making some of them massive. This effect 
appears as follows: on the field theory side, the multi-trace interactions violate the 
conservation of the energy-momentum tensor of each CFT, they retain, however, the 
conservation of a total stress tensor. On the dual gravity side, a linear combination 
of the k original gravitons remains massless while the rest of the gravitons obtain 
non-zero masses of order 1/N. One can think of this effect as a Higgs mechanism for 
gravity ||. 

Coupling more than one gravitons together, or giving the gravitons a mass, have been 
long standing theoretical problems @, [|, ^, (?], ||, [9], [l(| (see [11] for a recent review). 



The issue at stake is the possibility of coupling non-trivially massless gravitons, or 
giving them a mass without rending the theory UV sensitive. So far there are no-go 
theorems for non-trivial couplings among massless gravitons flO[| . On the other hand, 
most attempts to write an effective low-energy action for massive graviton theories 
are plagued by serious problems - ghosts, instabilities and strong coupling problems. 

This behavior is not only characteristic of theories where a graviton has a mass term 
in the Lagrangian, but also more exotic cases where the massive graviton can be a 
resonance. This happens when gravity is induced on branes, with the most celebrated 
example being given by the DGP model ||12|| . Indeed it was observed in ]l3| , [14| ] that 
this theory exhibits a similar non-decoupling behavior as the standard Pauli-Fierz 
theory. It was subsequently shown that the theory becomes strongly coupled at 



hierarchically low scales [15, Moreover, as has already been analyzed in detail in 



II |8| massive graviton theories are generically unstable. Instabilities similar to those 



of the Pauli-Fierz theory also occur in brane- induced gravity as reviewed in [17]. 

In this respect, having an example with a UV complete, non-perturbative formulation, 
where these issues can be addressed, is important not only from an academic desire 
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to know if such theories exist, but also for potential phenomenological applications. 
Infrared modifications of gravity of this type could be useful, for example, in the 
resolution of the cosmological constant. 

(ii) It is interesting to ask if and how standard properties of string theory and gravity (at 
the classical and quantum level) are modified in a multi-verse of interacting worlds. 
One can imagine new qualitative features on the string theory worldsheet and new 
dynamics in spacetime with potentially useful applications - for instance, potential 



applications in cosmology (c./. [18] for a search of cosmological solutions in a multi- 
gravity theory). 

The above AdS/CFT example offers a concrete arena to study the possibility and 
structure of such effects. In the context of the AdS/CFT correspondence, it has been 



argued [19] that multi-trace deformations on the boundary theory re-arrange the 
string perturbation theory in the bulk leading to a non-local string theory (NLST) 
with non-local dynamics on the worldsheet. One would like to know the precise 
rules of such dynamics. In particular, we would like to know whether such theories 
really define a new universality class of string theories. In this paper, we will have 
the chance to address some of these questions in a set of low-dimensional examples, 
where string theory is solvable. 

1.2 Deformed matrix model products and non-critical NLSTs 

The example of product CFTs and string theory on a union of AdS spacetimes can be 
generalized if we consider the product of more general quantum field theories (QFTs). As 
long as each of these QFTs has a string theory dual, the dual of the product theory will 
be a string theory on a union of spacetimes. 

One of the first and best understood examples of holographic duality in string theory is 
the duality between the double scaling limit of large N matrix models and two-dimensional 
quantum gravity coupled to conformal matter with central charge c < 1, i.e. c < 1 non- 



critical strings (for a review see [^, 21] and references therein). This is a good context 
for some of the above questions, because string theory in these examples is solvable at all 
orders in perturbation theory. Unfortunately, this case will not allow us to address the 
issues raised in point (i) above - the spacetime theory is a theory of scalar fields, hence we 
cannot arrange for a multi-gravity theory in this context. It allows us, however, to address 
some of the questions in point (ii) in a precise manner. Higher dimensional AdS/CFT 
examples that involve a multi-gravity theory and the related issues in point (i) will be 



discussed separately in a companion paper |2^] (see, however, subsection [L^ below and the 
discussion in section || for a summary of the main results). 

The precise setup we want to consider in this paper is as follows. The "boundary" 
theory involves a product of k large- AT matrix models Yli=i -Mi. The simultaneous standard 
double-scaling limit in each factor Mi gives the holographic description of a product of 
c < 1 non-critical string theories JX_i Si, where each factor in this product is independent 
and does not communicate with the rest. We want to deform the product rii=i-^i by 
adding multi-trace interaction terms to the total matrix model Lagrangian that couple 
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different A4j's together. The non-trivial question is: under what circumstances can we 
find new double-scaling limits with non-vanishing multi-trace interaction couplings? Such 
limits would define holographically a multi-verse of interacting c < 1 non-critical string 
theories. 

In section |2| we analyze the product of two Hermitian one-matrix models deformed by 
a general double-trace deformation. We will find that in this case a one-parameter family 
of double scaling limits exists with a non-trivial coupling between the two matrix models 
provided that the scaling properties of the single-trace operators that participate in the 
deformation are the same and that we tune the double-trace parameters to a special set of 
values. The free energy F of the deformed theory is no longer the direct sum of the free 
energies F\ and F2 of matrix models M\ and M2 respectively. It is the logarithm of a 
bilateral Laplace transform 

/oo 
dt + e M++-Fi(sin0 t++cos<9 t- )+F 2 (cos 9 t+-sin0 i_) ^ 
-oo 

where 9 G [0, 2tt) is a free angular variable parametrizing a £7(1) subspace of the three- 
dimensional double-trace coupling space and t+,t- are scaling parameters in the modified 
two-matrix model. This result is a natural generalization of a similar formula that arises in 
double-trace deformations of single Hermitian one-matrix models 24]. By holography, 
it gives the free energy of two coupled minimal string theories. 

In section ||| we find that the genus expansion of correlation functions in the deformed 



theory (1.1) boils down, when expressed in terms of correlation functions in the undeformed 
product M\ <S> M2, to a diagrammatic expansion where one has to sum over a series of 
terms that involve surfaces with contact interactions between the world-sheets of theory 
1 and 2. This makes the world-sheet theory non-local. At tree-level, where ( |1 . 1| ) reduces 
to a Legendre transform, all the contributing world-sheets are genus zero, however, at 
any higher loop order g > 1 there is a series of contributions from touching surfaces with 
different geni less or equal than g. This is a concrete example of a non-local string theory 
expansion along the lines of |19| . 

In higher dimensional AdS/CFT systems, in a limit where string theory in the bulk 
reduces to semi-classical gravity (this is the limit of large 't Hooft coupling for a gauge 
theory on the boundary), there is, at tree-level, an alternative reformulation of the gravity 
dual of multi-trace deformations as ordinary gravity with mixed boundary conditions for 



the dual fields g|, |§, ||, [28 



A similar reformulation of the tree-level theory appears to be possible in the context 
of matrix models and non-critical strings. This is in effect an old observation of Klebanov 
1 23], who pointed out that the string susceptibility exponents and tree- level correlation 



functions in the double-trace deformed theory can be interpreted in the dual string theory 
as a change of the Liouville dressing of the tachyon condensate from the right branch to 
the wrong branch. At face value, this is a peculiar statement in Liouville theory. It would 
seem to suggest that we have to change Liouville theory in a drastic way by modifying 
the boundary conditions of the tachyon field at infinity. This is at odds with our current 
understanding of Liouville theory, where both branches (or equivalently both the standard 
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and the dual cosmological constant interactions) have to be present in order to account for 



the right structure of correlation functions |3(j, [3l|, 33, 34 1. To the same effect, the only 
solution of the Wheeler-DeWitt equation that is regular in the strong coupling region has 
weak coupling asymptotics with both branches turned on pi]] . From this point of view, 
we can change the weak coupling asymptotics at the expense of introducing a singularity 
in the strong coupling region. The singularity could have a physical origin, e.g. it could 



be attributed to a bunch of ZZ branes [35] localized in the strong coupling region. It is 



unclear, however, why ZZ branes would suddenly make an appearance in this context. 
In section |] we will propose a different interpretation of the observations in [23], which 



agrees with the current understanding of Liouville theory and analogous observations in 
the AdS/CFT correspondence. The basic point is this. In Liouville theory the cosmolog- 
ical constant /i and the dual cosmological constant \x are both present and have a fixed 
relation. The theory is not modified if we exchange at the same time /i «-> jl and the 
branch of the Liouville interaction. However, each of these transformations separately will 
give a different theory. We propose that this theory is what the modified matrix model de- 
scribes at tree-level. In this modification Liouville theory continues to obey the usual rules. 
What we modify is the definition of the scaling parameter or in other words the specifics 
of the bulk/boundary dictionary. This is precisely what happens also in the AdS/CFT 
correspondence when we talk about mixed boundary conditions. 

Beyond tree-level the standard free energy evaluated as a function of fl does not agree 
with the exact results obtained from the Laplace transform fll-ip . This implies that the 
modified theory in the bulk is truly in a new universality class of string theories with a 



non-local worldsheet theory as envisioned in [19]. These points will be discussed further in 
section |] (and appendix ^]). 

There are several extensions of the modified product of two Hermitian one-matrix 
models that we explore in the main text (section ^ and appendices |A[ ^). One is the 
extension to matrix quantum mechanics and the dual bosonic c = 1 string theory. Another 
is the possibility to couple more than two non-critical string theories either by interactions 
that couple pairs of string theories or by higher-order 'vertices' that involve higher multi- 
trace interactions. As an illustration, appendix ^ considers double scaling limits in a 
product of three Hermitian one-matrix models coupled by a triple-trace interaction. 

1.3 AdS/CFT and multi-gravity 

A direct analog of a double-scaled matrix model in higher dimensions is a conformal field 
theory. In generic products of CFTs deformed by multi-trace interactions, the multi-trace 
couplings break conformal invariance. If $ = Y\j=i ®i ls the perturbing operator in terms of 
a set of single-trace operators Oi, where i labels the i-th CFT, then there is an upper limit 
on n in order for $ to be perturbatively relevant or marginal. This upper limit depends 
crucially on the spacetime dimension d through the unitarity bound. For scalar operators 
Oi the upper limit is 2 for d > 6, 3 for d = 5,4, 5 for d = 3 and infinity for d = 2. 

For concreteness, let us consider the case of two CFTs deformed by a double-trace 
interaction. Analyzing the one-loop (5- functions of single-trace and multi-trace couplings 



in conformal perturbation theory one recovers the following picture 22 
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At tree-level, i.e. to leading order in the 1/N expansion, there is a one-parameter family 
of non-trivial fixed points provided that the single-trace operators Oi participating in the 
double-trace deformation have the same scaling dimensions (or scaling dimensions with 
a sum equal to the spacetime dimension). This is one of the features that we encounter 
also in the corresponding matrix model analysis. In addition, one finds that the fixed 
points are repellors of the renormalization group (RG) equations, and perturbing away 
from them either drives the theory towards a fixed point, where the coupling between 
the CFTs vanishes, or towards strong coupling outside the range of validity of conformal 
perturbation theory. 

The dual description of this system involves quantum gravity on the union of two AdS 
spaces with mixed boundary conditions for the scalar fields that are dual to Oi jj], At 
tree-level, there are two massless, non-interacting gravitons in this system and the RG 
running of the double-trace couplings on the boundary is encoded subtly in the mixed 
boundary conditions (it is not visible, in particular, as a radial running of the background 
solution away from AdS). 

1/N effects modify this picture in an interesting way. On the gauge theory side, 
1 /N effects shift the submanifold of fixed points and induce an RG running of single-trace 
operators. On the gravity side, a non-trivial potential is generated for the gravitons leading 
to an interacting bi-gravity theory in the general spirit of Kogan and Damour . A linear 
combination of the gravitons remains massless, the orthogonal one obtains a mass of order 
1/N ||. The quantum generated potential backreacts to the original AdS solution and 
for generic bare values of the double-trace couplings the bulk solution is no longer the 
union of two AdS spaces. It is rather a union of two spaces that are radially deformed and 
interpolate between AdS spaces. 

In cases, where we can trust the boundary conformal perturbation theory, the gauge 
theory analysis predicts that the IR geometries are a union of AdS spaces with trivial 
coupling, hence in this region of spacetime there are again two massless, non- interacting 
gravitons. In a sense, the theory is dynamically washing away the mass of the graviton. The 
only way to retain a non-trivial bi-gravity theory everywhere in spacetime is to fine-tune 
the bare values of the double-trace couplings on a special one-dimensional submanifold of 
fixed points. The general lesson seems to be that massive gravity theories can in principle 
exist in a non-perturbative, well-defined quantum gravity context, but are not generic and 
require a high degree of fine-tuning. For a more detailed discussion of the general AdS / CFT 



case we refer the reader to ]22|. 

In the rest of this paper, we will deal with a solvable toy model of this general setup 
that involves matrix models and non-critical strings. The analogies between the matrix 
model and AdS / CFT setups will be summarized in the concluding section |5[ 

2. Matrix models and a multiverse of c < 1 string theories 

Multi-trace deformations in the context of matrix models and their implications for Liou- 



ville gravity were discussed originally in [37, 23, 24]. In this section, we will consider the 



effect of multi-trace deformations in a product of matrix models. Extending the analysis of 
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1 24 ] , we are looking for new double scaling limits that define holographically an interacting 
union of c < 1 non-critical string theories. As a first concrete illustration, we will discuss 
the case of two large- N Hermitian matrix models coupled by a double-trace deformation. 
Possible generalizations will be discussed in the last subsection ^3| and appendix [B]. 

2.1 Hermitian matrix models and minimal string theories 

Let us begin by recalling some basic facts about multi-critical Hermitian one-matrix models 
and their double scaling limits. The partition function of the k-th multicritical one-matrix 
model is defined by the matrix integral 

with 

k 

V k (<S>) = ^2(-y +1 c^ 21 (2.2) 
1=1 

and Cj a set of known constants [38]. In ( |2.1| ) we have chosen to consider a deformation 
of the potential by a term proportional to the single-trace operator Tr$ . This model 
has a double scaling limit, where TV" — > oo and ci + A — > with the scaling variable 
t ~ (c2+A)iV 2fc /( 2fc+1 ) kept fixed. In this limit, the partition function Tk = log becomes a 
function of the scaling variable t and the matrix model provides the holographic formulation 



of the minimal (2, 2k — 1) bosonic string [38]. 



We will consider a two-matrix model, which arises from the direct product of a k-th 
and a p-th multi-critical matrix model M.2,2k-i ® -M.2,2p-i by a double-trace deformation. 
The standard double scaling limit of this theory at zero double-trace coupling describes 
the decoupled union of a (2, 2k — 1) and a (2, 2p — 1) minimal string. 

2.1.1 A prototype: deforming the 7^2,3 ® A^2,3 product 

The theory we want to solve is a double-trace deformation of the A-^2,3 <S> M-2,3 two-matrix 
model with partition function 

Z = j ^^e^iTr^+Ai^j-^fi^+Aa^ _ 

(2.3) 

The double-trace deformation involves the operators Tr<£^, Tr^. It is g\2 that couples 
the two separate matrix models, but g\ \ and 522 will play an important role when we take 
double scaling limits. The ranks of the matrices $i, ( I > 2 are N\,Ni respectively. In this 
paper we use the notation 

Ni = N , N 2 = vN . (2.4) 

In the large iV limit both N\ and N2 will be scaled to infinity; the ratio v will be kept fixed 
and will be treated as an extra parameter of the system. The single-trace and double-trace 
coupling constants are such that the overall Lagrangian scales like N 2 . 
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It will be convenient to re-express the double-trace deformation as a sum of two squares 
<7ii (Tr $i) 2 + 5 22 (TV $l) 2 + 2512 Tr $f Tr $j = 

= n (cos6lTr^ + sin6lTr^) 2 + r2(-sin6lTr^ + cos6ITr^) 2 . (2.5) 
Setting for quick reference 

C = cos#, S = sinO (2.6) 
we deduce the double-trace coupling relations 

511 = nC 2 + r 2 S 2 , 522 = nS 2 + r 2 C 2 , 512 = (n - r 2 )SC . (2.7) 

With these definitions 

det 5 = 511522 - 5i2 = H72 , Tr 5 = 511 + 522 = n+r 2 . (2.8) 

We shall distinguish between two different cases: (i) T\T 2 7^ and (ii) r\r 2 = 0. In 
the first case, the double-trace deformation consists of two quadratic terms. In the second 
case, there is a single quadratic term or no deformation at all. 

• Case {%): r\r 2 7^ 

Let us define the free energy T of the single-trace theory as 

e H\N*] = j m e -NTr (2 Q) 

A common trick to deal with the double-trace deformations is to use the identity 



e 9U = -== / dy e ~~e nyu . (2.10) 
V4vr5 



In the case of (2.3), this trick allows us to recast the partition function as a double integral 
over single-trace parameters with a specific Gaussian weight 

Z= NlN2 r^dyae^^-^+^^+^^-^^+^^e^S+l] . (2.11) 
kKV^Jrxr 2 7_oo 

It is now possible to rewrite Z'm.a, more explicit form thanks to the well-known solution 
of the one-matrix model 38]. For the k = 2 multi-critical matrix model, in particular, 
the free energy reads 

JF[A, N 2 } = N 2 (-aix + \a 2 x 2 ^j + F(x, N 2 ) , (2.12) 

where 

F(x,N 2 ) = N 2 (~a 3 x 5 / 2 + ...) +^°(-^ log * + •••) + 

+N- 2 (a 4 x- 5/2 + ..) + 0{N~ 4 ) , (2.13) 
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with 

x = c 2 + A , 01 = 4 , a 2 = 576 , a 3 = 6144\/3 , c 2 = — , • • • (2.14) 

48 

The omitted terms inside each parenthesis in the expression of F ( |2.13 ) will be subleading 



in the double-scaling limit at each order in N and will not play a role in our discussion. 
Our case involves the free energy of two decoupled one-matrix models with parameters 



xi, x 2 . These parameters, which are to be integrated over in ( 2.11 ), are related to the 
variables y\ , y 2 by the linear transformation 

yi = Ai - Cx\ - Sx 2 , Ai = CAi + uSA 2 (2.15a) 

2/2 = A 2 + 5xi - Cx 2 , A 2 = -5Ai + ^CA 2 , (2.15b) 

where we define 

Ai = c 2 + Ai , A 2 = c 2 + A 2 . (2.16) 



Inserting this information into ( 2.11] ) and dropping an unimportant overall constant 



factor we obtain the partition function 

/OO AT 2 [Y „ A 1 C , A 2 S \ , ( n A 1 S A 2 C \ 1 

dx!dx 2 e-^- 2ai -^ + ^> 1+ (- 2aiv -^-^r) X2 \ x (2.17) 



-OO 

" 2 [(2a a +^+^)x?+(2o a +^+^)x3+2CS(i-i)xix a ] ^(xi.JV? )+F(x 2 / i /,AT|) _ 

To diagonalize the quadratic term in the exponent of the integrand we rotate from (x\, x 2 ) 
to (x+, xJ): 

x x = Ufx + + U{x- , x 2 = C/ 2 + x + + U^x- , (2.18) 

where 

Uf = S, t/f = C , ?7 2 + = C , CA7 = -5 . (2.19) 
In terms of the rotated variables 



r*00 

N 2 I dx + dx^ e ^(E+ X++ E-^) e -N* m %xl-N* m ^ 



X 

e F{Utx++U^x-,N , i)+F({U+x + +U^x.)/u,Nl) 

where 

_j_ / AiC A 2 S\ , / Ai5 A 2 C\ ■ 



(2.20) 



-E = -2ai — + — =- H7f - 2aii/ + — !— + — =— C/? , (2.21a) 

V n r 2 J V r i r 2 / 

m| = -Ifl + l+4a 2 ±^^') . (2.21b) 
8 \n r 2 nr 2 / 

Depending on the sign of the saddle point mass squared parameters we can define 
a variety of double scaling limits. We will distinguish between the following alternatives: 
(a) both m± are positive, (6) one mass 2 is positive, the other is zero, (c) both masses are 
zero, or (d) at least one mass is tachyonic. 
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(a) occurs if and only if 



11 , 
— < -2a 2 , — < -2a 2 . (2.22) 
ri r 2 



In this case, the only sensible double scaling limit that we can take requires 

E ± 

N -> oo , so that ^-N A / 5 a7 5 = t± is kept fixed . (2.23) 

4m± 



After standard manipulations (see |24}| ), we deduce the partition function 

z = e F(c/+t ++ c/ 1 -t_)+F( i ,-V 5([ /+ t++C /- t _)) = e F( tl )+F(t 2 ) = Z{ t x ) Z (t 2 ) (2.24) 

with the obvious definition of the scaling parameters t\ and i 2 . The resulting theory is 
the product of two undeformed, decoupled matrix models, which describe the union of two 
decoupled (2, 3) minimal strings. The double scaling limit has driven the theory back to 
the undeformed point where the effects of the double-trace deformation are washed away. 
A similar effect in the case of a single matrix model was observed in p4| . 

This behavior is reminiscent of what happens with an irrelevant perturbation in a 
higher dimensional quantum field theory when we follow the renormalization group flow 
towards the infrared. In our matrix model example, going towards the infrared is achieved 
by the double scaling limit which zooms around a critical point. When the condition 
Q2.22 ) holds, the effects of the double-trace deformation fl2.5| ) disappear in the double 



scaling limit and the perturbation behaves as an irrelevant operator. In fact, this is more 
than an analogy: in higher dimensions there is a range of parameters where double-trace 



perturbations generalizing ( |2.5D are indeed irrelevant in the usual RG sense [22]. 
Case (6) is more interesting. Now m+ = and m 2 _ > 0, which is equivalent to 

— + 2a 2 ] (— + 2a 2 ] = . (2.25) 
n ) \r 2 ) 

This condition is satisfied by a two-parameter family of deformations. Notice that under 
(2.25) it is impossible to achieve a non-vanishing g\ 2 without turning on at the same time 



the other two couplings gn,g 22 - 

To see what happens to the free energy, we set m? + = in ([2.20 ), and take the double 
scaling limit 

N -> 00 , E ± , t+ = x + N 4 / 5 al /5 fixed , 



"3 

t_ = J^N^aJ 5 fixed , t+ = ^-N 6 / 5 a~ 2/5 fixed . (2.26) 



E~ 2/5 r E^ 



After standard manipulations the double scaled partition sum becomes 

dt + e t+t + +F(U+t + +U^t-)+F(v-V5(U+t ++ U^-)) _ (2 _ 27) 

-00 

Z depends now on the double scaling parameters and parametrically on 9, which is 

one of the deformation parameters (|2.5|) . 
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Eq. ( |2.27 ) is an interesting exact formula that deserves a few comments. First, we 
observe that the old scaling parameter t + has been transmuted to a new scaling parameter 
t+ . t+ requires scaling with iV 4 / 5 , whereas t + requires scaling with jV 6 / 5 . A similar change 
of critical behavior occurs in a single hermitian matrix model deformed by a double-trace 
operator as was pointed out originally in p3| , |24| ]. In the single matrix model case with 
a double-trace deformation of the form g(Tr<I> 4 ) 2 , the new critical behavior occurs when 
g = — Tji- and exhibits the partition function 

/oo 
dt e tt+F{t) . (2.28) 
-oo 

Eq. ( 2.27| ) generalizes this effect to the two-matrix example (|2.3|). 

Both ( |2.27| ) and ( 2.28| ) are bilateral Laplace transforms of the partition sum of the 
original undeformed theory. For non-vanishing g\ 2 (2.27) is not anymore the partition 
function of the decoupled product of two one-matrix models. The dual minimal string 
interpretation of this formula will be discussed in section ||. For g\ 2 = we expect to 
recover the partition function of a decoupled product of theories. This partition function 
is factorizable with one factor being the double scaled partition function of an undeformed 
one-matrix model and the other being the Laplace transformed partition function as in 
( 2.28 ). Indeed, we can check that when gi2 = 

Z = e F{v-^) / dt+e t + t ++ F(t + ) _ (2 _ 29) 



oo 



The non-perturbative meaning of eqs. ( 2.27] ) - ( |2.29| ) is not completely clear. For 
instance, there are well-known problems in defining the one-matrix free energies F(t) non- 
perturbatively. Presumably, we can avoid this problem by looking at a different class 
of matrix models, which are non-perturbatively well-defined (see e.g. [40] and references 
therein). Even then, however, we have to check the convergence of the integral appearing 
in the Laplace transform or whatever generalizes the Laplace transform. Despite these 
important potential issues, eqs. (2.27) - ( 2.29| ) work well at any order in perturbation 
theory, which is what we will focus on in this paper. 

As a final comment on this case, we point out that there is an interesting analog 
of the critical behavior captured by eqs. fl2.27| ), (2.29) in higher dimensional AdS/CFT 



examples [22]. In conformal perturbation theory one finds again a one-parameter family of 
interacting fixed points with non-vanishing double-trace couplings, which translate in the 
AdS/CFT correspondence to a system of coupled string theories on a union of AdS spaces. 
When both masses in ( |2.21b| ) are taken to be zero (this is case (c) above) we obtain 

g 12 = and g n = g 22 = -— or . (2.30) 

2a 2 

With gn = g 22 = — ^ we are dealing with the decoupled product of two one-matrix 
models tuned to their individual double-trace deformed critical points that were analyzed 
in H (see eq. (jj| )). 

Finally, when at least one of the masses m± is tachyonic, the theory goes into a 
branched polymer phase where we cannot define anymore a double scaling limit describing 
a sum over continuous surfaces 
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• Case (ii): r\T2 = 

Similar manipulations can be performed when the determinant of the matrix of double-trace 
parameters gij vanishes, i.e. when r\ri = 0. In this case, there is only one integration over 
single-trace deformation parameters in ( |2.11|) and one can obtain again different critical 
behaviors depending on the sign of the saddle point mass squared m 2 . In the stability 



region m > 0, we recover (2.24), the decoupled product of two one- matrix models. When 
the mass m is zero, which occurs for a one-parameter family of double-trace deformations 
parametrized by the angle 8, we recover the interacting product of one-matrix models 



Q2.27 ). When the mass is tachyonic the theory goes into a branched polymer phase. 



2.1.2 Comments on the general M. 2 2k-l ® -^2,2p l product 

Much of what we said above about the A^2,3®^2,3 product goes through unchanged to the 
more general M.2,2k-\ ® -M-2,2p-i product of one-matrix models, although some important 
changes in the conclusion occur when k ^ p. In the general case, the double-trace deformed 
partition function reads 



Z = J D$ 1 D<I>2 e "^[ Tryfc($l)+{c2+Al)Tr *i +Tryp(<I>2)+{c2+A2)Tr<I, 2] x 

e - [gu(Tr $4)2 +922(Tr $4 )2+29i2 Tr $4 Tr |] 



with the potential V& still defined as in ( |2.2j ). 

We recall that the free energy of the k-th. multicritical matrix model is 

= N 2 f-aix + \a2X 2 ^j + F k (x, N 2 ) , (2.32) 

where as before x = C2 + A, but now more generally 

Ffc(x ' = " 2lTT a3iv2x(2fc+1)/fc + - (2 - 33) 

Because of the k dependence of the non-singular part of the free energy Fk(x, N 2 ), the 
double scaling limit in the single-trace theory is 

N -> oo , x -> , t = x N 2k l {2k+ V fixed . (2.34) 



On the other hand, since the singular part of the free energy in ( [2.32 ) is /c-independent, all 



the manipulations regarding this part go through unchanged as in the .M 2,3 <S> ^2,3 case. 
In particular, the criteria for stability or instability remain the same. 

Important differences occur when we take double scaling limits. There are no substan- 
tial changes when k = p, in which case the expressions we derived in the ^2,3 ® ^2,3 case 
remain true with the obvious modifications of double scaling limits. On the other hand, 
when k 7^ p we find that that there is no sensible double scaling limit that leads to the 
analog of eq. ( [2.27 )■ The reason for the absence of a sensible double scaling limit is the 



difference between the scaling behaviors ( 2.34| ) for k ^ p. We observe a similar effect in 



conformal perturbation theory in higher dimensional QFT examples [22]. New fixed points 
with non-zero double-trace coupling g±2 exist only when the single-trace operators 0i,02 
participating in the double-trace deformation have the same scaling dimension. 
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2.1.3 More general deformations of the .M2 2fc-i <£> -^2,2fe-i product 

For a single k-th matrix model new double scaling limits can be defined by adding double- 
trace deformations gO 2 with a general scaling operator O after the necessary fine-tuning 



of g [24]. In that case, one obtains the free energy 

/oo 
dto e toto+Ht ' to) . (2.35) 
-oo 

This example can be generalized to 

n ;>oo 

f({t},{T}) =log[| / dti eE7=x*i*i+^({*}.OT) ; (2.36) 

where {T} is the set of coupling constants that remain unintegrated. 

Accordingly, for a product of A;-th matrix models new double scaling limits can be 
defined by tuning a set of multi-trace parameters. Consider, for instance, two fc-th one- 
matrix models deformed by the double trace operator g\\0\ + 522C2 + 2912O1O2, where 
0i, O2 are respectively scaling operators in the matrix models 1 and 2. We expect new 
double scaling limits leading to the free energies 

/oo _ , 

dT + e T+T++f k (tl,U+T + + U 1 r_)+^(t 2 ,C/ 2 r + +U 2 r_) ^ ^37) 
-00 

where t\,t2 are single-trace couplings for the lowest dimension operators. Further general- 
izations can be envisioned by introducing more couplings. 

2.2 Matrix quantum mechanics and c = 1 string theories 

Analogous statements can be made for a pair of matrix quantum mechanics theories coupled 
by a double-trace deformation. In appendix |A| we consider in some detail double-trace 
deformations involving the cubic single-trace operator Tr<3? 3 . The partition function of the 
theory we analyze there is 

Z = J D*i(t)£>$ 2 W e - 7V /o^ rf *[ Tr (l*?+l*?- A i*?)+ Tr (5*i+5*2- A 2*|)] x 

e ~ ft R ^[9ii(Tr1>?) 2 +2s 12 Tr$?Tr$|+<; 22 (Tr$3) 2 ] ^ 

The matrix model lives on a compact one-dimensional spacetime with radius R. Besides 
this extra feature (that can be dealt with a Fourier transformation) most of the elements of 
the analysis of the previous subsections go through also in this case. The possible double 
scaling limits and the resulting expressions are analogous to the zero-dimensional matrix 
model case above, so we will omit a detailed discussion here. The interested reader can 
find a detailed analysis of the partition function (|2.3§| ) in appendix |A[ 

2.3 Triple intersections and beyond 

So far we have restricted our attention to pairs of matrix models. In a suitable double 
scaling limit, these models provide the holographic description of a product string theory 
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with two components. There is a subtle communication between the individual components 
of this product induced by the double-trace deformation of the matrix model dual. There 
are various possible generalizations of this picture. 

We could consider, for instance, n one-matrix models coupled together by a general 
multi-trace interaction. We can imagine a large variety of possibilities. One of them is to 
take n theories coupled two-by-two via double-trace deformations. The overall theory is 
an n-matrix model with an action of the form 

n n 

S ^.S', • £ mp.Oj • (2.39) 

where Si are the single-trace actions and Oi are single-trace operators, e.g. Oi = Tr<I>|. 
Recasting gu in terms of n eigenvalues and an 0(n) matrix parametrized by n orthogonal 
unit vectors v a , we can write S as 

n n / n \ 2 

5 =E 5 *+E r « E^°0 • (2 - 40) 

i=l a=l \i=l / 

The corresponding partition function can be analyzed as above and allows for an obvious 
extension of the double scaling limits presented in previous subsections provided that the 
operators Oi have compatible scaling behaviors. 

Other possibilities include higher multi-trace interactions with three or more single- 
trace components. This allows for a more intricate set of double scaling limits. As an 
illustration, in appendix ^ we analyze the triple intersection of three 2nd multicritical 
matrix models with the action 

n -.3 

S = Y, S * + n ^ ^(Tr^)(Tr^ )(T rcI>4). (2.41) 

i=l i,j,k=l 

For example, in the special case, where 5123 is the only non-vanishing triple-trace coupling, 
we can find a double scaling limit with partition function 

/oo 
dhdt 2 g^i+^+ELi ^1+^+^*3) . (2 .4 2) 
-00 

U\ are constants that can be determined. Once again we obtain a modified partition 
function related to the undeformed one by a Laplace transform. A more intricate pattern 
of double scaling limits is expected for quartic and higher multi-trace deformations. 

One qualitatively new feature in the triple-trace example, compared to the double-trace 
case, is the possibility to find a non-trivial fixed point with 5123 ^ and no other multi- 
trace couplings turned on. This was not possible in the case of double-trace deformations, 
where at points with gn = 522 = 0, g\2 / at least one of the saddle point masses m± was 
tachyonic. 

In higher dimensional QFTs unitarity and the requirement that the n-trace interaction 
is a relevant or marginal operator puts an upper bound on n. For example, in d = 4 
dimensions the bound is n = 3. The bound is absent only in two or lower dimensions. The 
general structure of the one-loop /^-functions for multi-trace interactions in d dimensions 
will be discussed in [22]. 
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3. Correlation functions in coupled minimal string theories 



In this section we will discuss scattering amplitudes in a theory of coupled minimal strings 
using the matrix model definition fl2.27|) . We will illustrate the main points by recasting 
the correlation functions of the deformed product theory as a diagrammatic expansion in 
terms of correlation functions in the undeformed theory. As an example, we will consider 
two-point functions at tree-level and one-loop. The analysis in this section is a simple 
extension of the corresponding study of correlation functions in a single modified minimal 
string in . It will provide a concrete illustration of the non-local string theory structure 
underlying the dynamics of the string theories of interest and will be a good guide for the 
higher dimensional AdS/CFT examples discussed in the introduction, where it is much 
harder to explore this structure explicitly beyond tree-level. 

To be definite, let us focus on one of the cases discussed in section ||: two (2, 2k — 1) 
minimal strings coupled on the matrix model side by a double-trace deformation. The 
partition function of the modified theory is given by the Laplace transform 

/oo 
^ e \- 1 i + t + +F 1 (St + +Ct^,t 1 )+F 2 {Ct + -St-,t 2 ) _ (3 ]\ 

-oo 

t\,ti are single-trace couplings for the scaling operators 0\, O2 of theory 1 and 2 respec- 
tively, whose correlations functions we will compute. These operators do not partake in 
the double-trace deformation that couples the two minimal strings. The constant A has 
been inserted in ( |3.1| ) to help us keep track of the genus expansion of the free energies. 

We can obtain the all-genus diagrammatic expansion of correlation functions in the 
deformed theory by expanding the integral expression ( |3.l|) around the saddle point value 
of t+, which we will call t s . t s is determined by the following equation 

t+ + S{Vi)\ st a +ct- + C (7 3 2)| ct a -St- = , (3.2) 

where the subindex next to the one-point functions denotes the point at which the corre- 
lation functions are evaluated. In what follows, we will keep this index implicit. V\ and 
V2 are operators in theory 1 and 2 respectively, whose single-trace couplings St + + Ci_, 



Ct+ — St- appear inside the integral in ( |3.1| ). In (|3.2| ) we used the n = 1 version of the 
identities 

\- 1 {V?)\ t = d?F i {t), i = l,2. (3.3) 



The one-point functions in ( |3.2[) receive contributions at all geni. Since we are inter- 
ested in a genus expansion (in other words, an expansion in powers of A) of F, it will be 
useful to define our point of expansion by the tree-level version of ( |3.2; ) 

t + + S(P l )o + C(P 2 )o = . (3.4) 
The subindex in correlation functions will denote from now on the genus. 
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Setting t+ = t s + 1 for the integration variable in ( |3.l|) , we expand around t s to obtain 
the free energy expression 

F{t + ,t-;ti,t 2 ) = X'H+U + F^Sts + Ct-,h) +F 2 (Ct s - St_,i 2 ) + (3.5) 



+ log / dt exp 



oo oo oo 

t £ a 9 - 1 {S{v x ) g + c(v 2 ) 9 ) + E E i^ 9-1 ( sn ( v i) 9 + cn ( p %) 9 

g=l n=2 g=0 

For example, from this expression we read off the tree-level and one- loop free energies 

Ft r ee(h,t--,h,t2) = X-%t s + Fi 0) (St s + Ct-,t 1 )+F^\ct s -St-,t 2 ) , (3.6a) 

F one _ loop (t + ,t.;t 1 ,t 2 ) = F 1 (1) (St s +Ct_^ 1 )+F 2 (1) (Ct s -St_^ 2 )-ilog [5 2 (PiPi) + C 2 (V 2 V 2 ) } 

(3.6b) 

where Fi = Yl'gLo X 9 ~ 1 F^ S (i = 1,2) is the genus expansion of the free energy of each 
theory. We see that the tree-level expression ( |3.6a| ) is merely a Legendre transform of the 
total undeformed free energy. On the other hand, the one- loop expression ( |3.6b| ) receives 
contributions both from tree-level and one-loop quantities in the undeformed theory. 

At any genus, correlation functions can be computed with the use of the Feynman 
diagrams in fig. [l|, which include vertices (tadpoles, 2- and higher n-point vertices) of both 
minimal strings. Vertices of either string can be connected with a common propagator. The 
strength of the vertices and the actual form of the propagator is controlled by the double- 
trace parameter 9. The generic diagram of the modified theory is a sum of correlation 
functions in the original product theory on disconnected worldsheets. One could reformu- 
late the new interactions in terms of a set of non-local interactions on the worldsheet of 
the original theory. Because of its tractability, our setup provides a concrete example of 



the non-local string theory construction proposed in [19|. We will return to this point in 
the next section. 

To illustrate the generic structure of correlation functions in the modified product 
theory with a few examples, we will now consider the two-point functions of the operators 
Vi, Oi (i = 1,2) at tree-level and one-loop. 

3.1 Two-point functions at tree-level 

At tree- level the free energy of the modified theory is the Legendre transform of the original 
theory. We can obtain the tree-level amplitudes of the scaling operators 0%, 2 by differ- 
entiating ( 3.6a| ) with respect to t\, t 2 and using the chain rule. For two-point functions we 



find (we will denote the modified 2-point functions as ((•••))) 

«<W>„ - <^>o + s2{v X^ V chv 2 v 2)o ■ ^ 

((o 2 o 2))o = { o 2 o 2)o + s2{v ^ITchv 2 v 2 ) • < 3 - 7b ) 

{{Ol ° 2))0 ~ *Wi>o + CW 2 >o • (3 ' 7C) 
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X 



8 




Figure 1: Vertices necessary for the diagrammatic expansion of the modified free energy of two 
coupled minimal strings. The white (grey) worldsheets are respectively worldsheets of minimal string 
1 (2). The vertices in the last line apply to any genus g = 0,1,2,... and can have n — 3,4,5,... 
external legs. 

As expected, the last correlation function that describes the non-trivial interaction between 
the two minimal strings is proportional to SC, which vanishes only when the double-trace 
coupling g\2 is zero. Diagrammatically, these two-point functions can be recast as 



«0101» O = Q + (7y^7^ , (3.8a) 



«020 2 »o = Q + 0~~0 ' (3 ' 8b) 

((O 1 O 2 )) = Q—Q . (3.8c) 

Differentiating the free energy F with respect to the scaling parameters i+, t- gives 
correlation functions for the operators V+^V-, which are defined in the following way. In 
the original theory, before the double-trace deformation, V± are operators in the product 
theory defined by the linear combinations 

r + = svi + cv 2 , V- = cvx - sv 2 . (3.9) 

The couplings of V± are what we called above t±. In the deformed theory (jjQ]), t + is 
modified to t + and i_ remains invariant. Accordingly, V+ is modified to V+ and V- 
remains invariant. 

At tree-level, V+ and V- exhibit the one-point functions 

((P+))o = ^ , ({V-)) = A- 1 (P-) (3.10) 
and the two-point functions 

"^''•^sWftTcw*' (3 ' lla) 

«*.*-»„ = cow* + s*{vw« - %^lc^! • P-ub) 



CS{(V 1 V 1 )q-(V 2 V2)o) 
S 2 (V 1 V 1 ) + C 2 (V 2 V 2 ) ' 



((^-))o = - i^^™ . (3.11c) 



3.2 Two-point functions at one-loop 

One-loop correlation functions can be obtained from ( |3.6b ) in a similar fashion. Here we 



will list the two-point functions of the operators Oi. Analogous expressions can be deduced 
for the amplitudes of the operators V + , V- ■ It will be convenient to set 

K = S 2 (VxVx)o + C 2 (V 2 V 2 ) . (3.12) 

The two-point function ((dOi)) 1 reads 

((o 1 o 1 )) 1 = (OiOi)i + s 2 ic- 1 (r 1 o 1 ) (r 1 o 1 ) 1 - s^^xiVxOx^ivfh - 

-C 3 5 r3 /C- 3 <Pi)i(^iOi)g(^) - CS 5 Kr*(V 2 )x(VxOx)l(Vl) Q - C 3 S 2 fC- 3 {V 2 )x{VxOx) 2 {Vi) + 
+S 4 K- 2 {V 2 l ) 1 {V l 1 )l + CS*Kr 2 (Vl)x(VxOx)l + S 2 fC~ l {Vx)x(VxO 2 ) + C5/C- 1 (P 2 ) 1 (P 1 O 2 ) + 

+is 4 £- 2 <p?0 1 ) o + 5 6 /C- 3 (P 1 2 1 )o(^iOi)o(Pi 3 )o + C7 3 5 3 /C- 3 (P 1 2 1 )o(PiOi)o(P 2 3 )o - 

As^iVfoDo - 5 4 /C^ 2 (P 1 1 ) (P 1 3 Oi)o - ^ 4 /C- 2 (P 3 ) (Pi0 2 )o - 

-ic5 3 /C- 2 (P|) (Pi0 2 )o - ^5 6 /C- 3 (PiOi)^(Pi 4 )o - \c A S 2 K,- z {VxOx)l{Vt)o + 
+5 8 /C- 4 (P 1 1 ) 2 (P 3 ) 2 + 2C 3 5 5 /C- 4 (P 1 1 )g(^i 3 )o(^!)o + C & S 2 Kr\VxO x )l{Vl)l (3.13) 
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♦ GP + G^O + CMZX> + O-Cp 
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Figure 2: TTie diagrammatic expansion of the one-loop amplitude {{0\0\)) l . 

and has the diagrammatic expansion in fig. |2[ A similar expansion, with obvious modifi- 
cations, applies to ((0 2 2 )} 1 . 

The expansion of ((0\0 2 )) l is 

((0i0 2 ))i = CSJC- 1 (V 1 O 1 ) 1 (V 2 O 2 ) + CSK- 1 {V 1 1 ) Q {V20 2 ) l - 
-C5 5 /C- 3 (P 1 ) 1 (P 1 1 ) (P 2 2 )o(^i 3 )o - C 2 5 4 /C- 3 (P 2 ) 1 (P 1 1 ) (P 2 2 )o(Pi 3 )o - 
-C 4 5 2 /C- 3 (P 1 ) 1 (PiOi) (P 2 2 )o(^!)o - C 5 5/C- 3 (P 2 ) 1 (P 1 1 ) (P 2 2 ) (p|)o - 
-C 2 5 2 /C- 2 (P 1 ) 1 (P 1 1 ) (P|0 2 )o - C 3 5/C- 2 (P 2 ) 1 (P 1 1 )o(P|0 2 )o + 
+C5 3 /C- 2 (P 1 1 ) (P 2 2 )o(^i 2 )i + C 3 5/C- 2 (P 1 1 ) (P 2 2 )o(P|)i + 
+C5 3 /C- 2 (Pi) 1 (P 2 Oi) (P 2 2 )o + C 2 S 2 JC- 2 (V2}i(V!0 1 }o(V20 2 }o + 

+\c 2 s 2 K- 2 {v 2 l o l ) Q {vlo 2 ) Q - ^cs 3 jc- 2 (v!o 1 ) (v 2 o 2 )o - ^C 3 5/C- 2 (P 1 1 )o(P|0 2 )o+ 
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+CS 1 K-\V 1 1 )v{V 2 2 UVl)o + C 7 5/C- 4 (P 1 1 )o(P 2 2 ) (P|)o + 



+2C 4 5 4 /C- 4 (P 1 1 )o(P 2 2 )o(^i 3 >o(^2 3 >o + C 2 5 4 /C- 3 (P 1 1 )o(P|0 2 )o(P 1 3 ^ 



-^CS 5 lC- 3 (V 1 1 )o(V202}o(Vl}o - ^C 5 SfC- 3 (Vi0 1 } (V20 2 )o(Vl)o . (3.14) 

The corresponding diagrams appear in fig. ||. 

Most of the terms appearing in equations (3.13) and ( |3,14| ) have a unique representation 
as diagrams in figures |2| and ||. There are, however, some exceptions. These include the 
last three diagrams in fig. Q and their analogs in fig. ^. For example, both diagrams 
in fig. U can be associated to the amplitude contribution (ViO^Q^Pf Viewing the 
contact interactions as tiny wormholes that connect different parts of a worldsheet, we 
recognize these diagrams as representations of the same amplitude in a different region of 
the worldsheet moduli. In string perturbation theory one sums over all these configurations 
automatically. 

There are many contributions to the one-loop renormalization of the two-point func- 
tions ( ggp - dOH )- 

Some of them are standard torus amplitudes in theory 1 or 2 sep- 
arately, but the majority comes from contact interactions where worldsheets of theory 1 
and/or 2 connect via a tiny 'neck' (wormhole) associated to a common propagator /C _1 
J3.12| ). Contact interactions between worldsheets of different theories are especially in- 
teresting since they provide a clear picture of how theories 1 and 2 communicate with 
each other on the level of the worldsheet. Such interactions contribute both to single 
theory two-point functions (((0\0\)) or ((0 2 2 ))) and to mixed two-point functions of 
the form {{0\0 2 )). In the latter case the 1-2 interactions work to renormalize the tree- 
level result ( |3.8d ) either by renormalizing separately the single theory two-point functions 
{0\Pi)o, (0 2 V 2 )o or by renormalizing the propagator /C _1 with a double wormhole inter- 
action C 2 S 2 Kr 2 {V\0\)^Vl0 2 ) (the last diagram in the third line of fig. |). A similar 
diagram renormalizes the mixed two-point function {g\g 2 ) for the two AdS gravitons gi,g 2 
in the AdS/CFT example of [jl|, [2j. In that case, this is precisely the effect that lifts the 
mass of a certain linear combination of the gravitons and leads to massive gravity. 



4. Worldsheet and spacetime aspects of the multi-trace coupling 

In this section we want to discuss the worldsheet and spacetime interpretation of the 
theories defined above through the matrix model. There are two aspects we want to 
emphasize: (a) the modified genus expansion and non-locality on the worldsheet, and 
(6) the analogy with higher dimensional AdS/CFT examples and the possibility of an 
alternative formulation of the deformed theories that does not involve non-locality on the 
worldsheet. 

4.1 Touching surfaces in Liouville gravity and the NLST structure 

In conventional matrix models, i.e. matrix models with only single-trace interactions, the 
Feynman diagram expansion generates discretized random surfaces from which the contin- 
uous string worldsheet arises in a suitable double scaling limit. Multi-trace interactions 
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Figure 3: The diagrammatic expansion of the one-loop amplitude ((OiC?2))i- 
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Figure 4: Two diagrams corresponding to the same amplitude contribution (VxO\) q(Pi)o 



modify this expansion by gluing together different random surfaces at a plaquette. In the 
continuum limit, one can think of this contact interaction as a tiny neck (wormhole) that 
creates a network of touching surfaces. Such microscopic degenerations of the string world- 
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sheet are already present in the conventional theory without the multi-trace deformation, 
however, they give a sizable contribution that modifies the conventional result only for spe- 
cial values of the multi-trace couplings where a new double scaling limit is possible. The 
diagrammatic expansion of the previous section gives a clear demonstration of this effect 
of contact interactions between different worldsheets. In fact, the examples of section || 
give an even more dramatic illustration of this general phenomenon: even worldsheets of 
different string theories can touch via a common wormhole. 

One can argue on general grounds that multi-trace interactions will have a similar effect 
also in the AdS/CFT correspondence in higher dimensions. String theories resulting from 



such a deformation were termed non-local string theories in [19], since one can presumably 
reproduce the effect of the tiny wormholes with a non-local deformation of the worldsheet 
action. More specifically, if N w is the number of disconnected components of the worldsheet 
s = eW e e (2) e • • • e £ (iVw) one can envision a non-local worldsheet action of the form 

S WS = Y,\ d 2 aWy/^)£ J d 2 af ) d 2 4 3) Vg((ai)g(<J 2 )g[X(a 1 ),X(a 2 )} + 

i=l *,i=i 
+trilocal and higher — order interactions (4-1) 

where X denotes collectively the target space fields and Q is a non-local interaction. Q 
and the rest of the higher order interactions in Q4.1| ) should be determined directly in 
string theory by consistency, e.g. by requiring the cancellation of the Weyl anomaly. The 
precise rules for such theories have not been worked out, however, the consistency of the 
AdS/CFT correspondence implies that such non-local string theories should exist. In fact, 
we recognize a concrete example of such a theory in the non-critical NLSTs defined in 
section ^. In that case, see section [| we can determine explicitly the precise rules for 
the diagrammatic expansion of the correlation functions and from them one could work 
backwards to deduce the required form of Q etc. 

4.2 An alternative interpretation? 

When a CFT admits a dual supergravity description, it has been argued [25, 26, 27, 28, 29] 



that there is an alternative way to understand multi-trace deformations as mixed boundary 
conditions for the dual fields in AdS. In the matrix model case, there are indications for a 
similar reformulation of the deformation on the string theory side as a local effect on the 
worldsheet that involves changing the branch of the Liouville dressing of a vertex operator 



1 23]. Here we will argue that there is a simple interpretation of this observation at tree- 
level, which is in accordance with our current understanding of Liouville theory and the 
analogous observations in the AdS/CFT correspondence [25, [2(1, 27, 25, EjJ. Beyond tree- 



level, a similar reformulation does not go through, and one has to think of the deformed 



theory as an NLST along the lines of [19]. 



4.2.1 The tree-level theory 

To illustrate the main point, let us consider first the simpler situation of a single one-matrix 
model deformed by a double-trace deformation. The analogous AdS/CFT example involves 
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a CFT deformed by a double-trace interaction of the form 

5S= [ d d x gO 2 , (4.2) 



where O is a single-trace operator with scaling dimension A. g is a coupling that scales 
like N°. It will be useful to recall first the main features of this example. 

If the scaling dimension A is less than |, the perturbing operator ( |4.2p is relevant 
and the theory runs towards an infrared (IR) fixed point, where O has a different scaling 
dimension d— A. At large N, the IR theory is simply related to the UV theory by Legendre 
transform |ft2]| . At tree-level and within the supergravity description, this RG running on 
the boundary has a holographic counterpart in (d + l)-dimensional AdS space as a flow 
between different boundary conditions for the bulk field <fi which is dual to the operator O. 

There are two points we want to emphasize with respect to this statement. First, in 
contrast with the case of a single-trace perturbation, where the dual gravity background 
is deformed away from AdS, the multi-trace deformation does not backreact to the bulk 
geometry to leading order in 1/JV. Second, changing the boundary conditions of the bulk 
field does not mean that we pick a different solution of the bulk equations of motion for 
the dual field <j>. This is still controlled by regularity in the bulk. Instead, we modify the 
definition of the source - in other words, we modify the bulk/boundary dictionary. This 
last point will be important for the matrix model discussion below, so we take a moment 
to explain it here in slightly more detail (for a more detailed exposition of this well-known 



material we refer the reader to 43], which is what we will follow here mostly, and the 
standard references therein). 

In coordinates where the Euclidean AdS metric takes the form 

ds dA 

consider a scalar field A with the action 



- l = ^(dr* + jydx i f\ (4.3) 



S = j d d xdr^g- Q W) 2 + ^V 2 j + J d d Xy/g C boundary \ r=e . (4.4) 

4> is the scalar field dual to the single-trace operator O. The mass m is related to the 
scaling dimension A via the relation 



A=^-u, u = ^ + m 2 . (4.5) 



We have included a boundary term to the action (4.4), which is defined at r = e, the IR 
bulk regulator. To reproduce the double-trace deformation (|4.2|) we are instructed to use 
the boundary interaction 

{--boundary — ^ ' (^'^) 



/ has a simple relation with the field theory g in (4.2). A careful treatment gives 



f = -A-2g^(2n^^^l) . (4.7) 
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Varying the action ( [4.4|) with respect to (f> we deduce the mixed Neumann/Dirichlet bound- 
ary conditions 

/ (f)(x, e) + d<j>(x, e) • h = , (4.8) 

where h = er is a unit vector specifying the normal to the boundary. 

To compute correlation functions in the AdS/CFT correspondence, we evaluate the 
on-shell bulk action as a functional of the boundary source (ftb- With mixed boundary 
conditions (f4,8j), the boundary value problem is 

(a-m 2 )(j)(x,r) = , (4.9a) 
/ </>(x, e) + d(j)(x, e) • n = <t> h {x) . (4.9b) 

After the Fourier transformation 

<f>(xi,r) = j^j^ j d d k e^Vte.r) (4.10) 

the wave equation (4.9a) becomes (set 4>(k,r) = r d / 2 x(k,r)) 

r-^j + fcV + v 2 1 x(k, r) = . (4.11) 

The unique solution to (PI) , (|4.9b| ), which is regular at r — > oo, is given by the modified 
Bessel function of the 2nd kind K, u . Consequently, 

r ) = ( f m vr^m \ ^ r > • M k > r ) = ^^(^) • ( 4 - 12 ) 

In particular, we notice that for any bare value of the double-trace coupling g, the solution 
is always the same Bessel function. What changes is the relation between the source fa 
and the asymptotic coefficients a, /3 in the near boundary (r — ► 0) expansion of (j) 

cj)(x,r)~r d - A [a(x)+0(r 2 )] + r A [j3{x) + 0{r 2 )] . (4.13) 

Regularity, or equivalently the fact that we chose IC U as the solution of ( [4.11 ), has already 
fixed a linear relation between a and (5. 

Having said this, let us return to the matrix model case and the dual minimal strings. 
The analog of the wave equation ( 4.11 ) is the Wheeler-DeWitt (WdW) equation. For a 
minimal string with Liouville interaction 

tS IM ->ff,a Mr *>», a+ =6>0 (414) 

the wavefunction ^f u for a mode that corresponds to the vertex operator {Q is the linear 
dilaton slope) 

V u = O matter e-^ +Q ^ = V v e-%* (4.15) 
obeys, in the mini-superspace approximation, the WdW equation 

l^\ 2 + 4fx£ 2 + u 2 U^) = 0. (4.16) 
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i is related to the zero-mode of the Liouville coordinate <po by the relation 

i = ei 6 ^ . (4.17) 

The weak coupling end of Liouville theory lies at I — * 0, or equivalently (p — > — oo. 

The similarity between eqs. (4.11) and ( 4.16j ) is obvious. Again, the only solution that 



is regular in the "IR", i.e. the strong coupling region at I — * oo, is 

oc K v {2^Jii) . (4.18) 
Writing K, u in terms of the modified Bessel functions of the 1st kind X v 

M -v-Tu) (4-19) 

we observe the weak coupling asymptotics 

K. V {1) ~ / ( t^—X + ... - V^—^ + ■ (4-20) 
2sm{irv) \L (1 — v) 1(1 + ^) / 

Hence, regularity of the full wavefunction requires the presence of both the 'right' branch 
with v < 0, that satisfies the Seiberg bound |44}| , and the 'wrong' branch with v > 0, with 
a fixed reflection relation. 

Deforming the matrix model with a double-trace operator, e.g. the operator (Tr<I> 4 ) 2 , 
one finds two possible double-scaling limits: the standard one at double scaling parameter 
g = 0, and another one with different string susceptibility exponents at a special non- 



vanishing value of g [23, p4]. These two critical points are analogous to the UV and IR 



fixed points of the CFT deformation in (|4.2|), The analysis of the string susceptibility 



exponents [23] and tree- level correlation functions |41] in the new double scaling limit 
indicates that one can reproduce the matrix model results in a Goulian-Li approach [45] to 
the minimal string by changing the Liouville dressing, from the right branch to the wrong 
branch, of the vertex operator that is dual to the matrix model operator involved in the 
double-trace deformation. Ref. proposed this as a possible interpretation of the string 
theory dual to the deformed matrix model. 

The statement about changing the branch of the Liouville dressing cannot mean that 
we modify the boundary conditions of the WdW equation in such a way that the new 
solution is solely 1\ v \- This solution is singular in the strong coupling region (£ — > oo), and 
one would need additional ingredients to explain this singularity, e.g. the singularity could 
be accounted for by the presence of a large number of ZZ branes |35f] localized at the strong 



coupling region (see [46] for a related discussion). It is unclear, however, why ZZ branes 
would all of a sudden appear in the minimal string as we change g continuously from the 
undeformed point to the new critical point. Also, this interpretation would not fit with our 
understanding of the analogous situations in the AdS/CFT correspondence, as described 
above. 

Along the same lines, in recent years it has been understood that in Liouville theory 
both the standard Liouville interaction ( 4.14| ) and its dual 



5S Liouv iii e = fx J d 2 z O mat tere a ~' p , a- = -^Q + ^j (4.21) 
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have to be present with a fixed relation /2 = /i(/x). Otherwise one cannot explain the 
structure of the exact two- and three-point functions |30], [3l], 32, 33, 34]. 

In view of these observations, it is more appropriate to think of the deformed theory 
at tree-level in the following way. As a specific example, let us consider the case of the 
2nd multi-critical matrix model deformed by g(Tr <1> 4 ) 2 at the new critical point where the 
partition function is given by ( 2.28| ). At tree-level, we can still think of the string theory 



dual of this matrix model as the usual (2,3) minimal string - the same as that at the 
undeformed g = point. We should, however, modify the bulk/boundary dictionary and 
think of the amplitudes of the new theory as a function of the dual cosmological constant fl. 
A way to rephrase this proposal is as follows. It is known that Liouville theory is symmetric 
under the simultaneous exchange of fj, <-> jl and a+ «-> a_ in ( [4.14 ), (4.21). Performing 



half of this transformation, e.g. /i <-> fi with a+ unchanged, gives a different theory. This 
is the theory we propose as the dual of the deformed matrix model at tree-level. This 
interpretation is consistent with our current knowledge of Liouville theory, the checks of 
the string susceptibility exponents and tree-level correlation functions in |23|, [24], |4l| and 
meshes nicely with the more general picture of double-trace deformations in the AdS/CFT 
correspondence. 

We propose a similar interpretation of the double scaling limit ( |2.26| ) that leads to 
the modified partition function (|2.27f ) of two coupled minimal strings. In that case, we 
re-interpret the original product partition sum as a function of the cosmological constant 
t- = Ct\ — Sti and the dual t + = St\ + Ct^- Some evidence for this proposal is provided 
in appendix |C] where the sphere partition sum is analyzed from the continuum formalism 
viewpoint. 

4.2.2 Beyond tree-level 

The full partition function of the theories we considered in this paper is given by the 
Laplace transform of the original partition function. For instance, the free energy of the 
standard double scaling limit of the 2nd multi-critical matrix model is 

F(t) = --t b/2 - — log t + — r 5/2 + 0(t~ 5 ) . (4.22) 
w 5 24 & 2160 y J v ; 

The Laplace transform of the exponential of this expression gives the free energy of the 
theory modified by (Tr<£ 4 ) 2 at the new critical point 

F (t) = |?/3 _ L Jogt + ZLt-5/3 + O(t ~-10/ 3) (4 23) 

5 3d 9dU 

The relation between the cosmological constant t and its dual t is, up to a constant 
numerical factor, 

t = P /3 . (4.24) 

This relation reproduces the scaling of all terms in the perturbative expansion, and can 
even be used to reproduce exactly the tree-level result with appropriate fixing of the nor- 
malization of t, as explained in the previous subsection. It fails, however, to reproduce the 
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higher loop coefficients. For example, it fails to reproduce the torus coefficient — ^ giving 
instead the factor — 

This implies that the tree-level interpretation of the previous subsection does not 
extend to higher loops, where one should think of the modified theory as a genuine NLST. 
It is tempting to speculate that a similar result holds in the higher dimensional AdS/CFT 
example, and that also there it is necessary to invoke the NLST structure to explain the 
properties of the modified string theory in the bulk. 

4.3 Comments on the spacetime properties of the modified theories 

The overall picture for the string theory dual of two double-trace deformed large iV matrix 
models at the critical point (2.26), ([2.27 ) is a pair of minimal string theories with a subtle 



interaction induced on the level of the worldsheet via tiny wormholes. It is interesting to 
ask how such interactions manifest themselves in spacetime, or in other words how they 
enter in the string field theory action. 

It is apparent from the scattering amplitudes computed in section ||| that there will be 
several types of contributions to the overall effective string field theory action, which will 
look schematically as 

Stotal = «Si + S 2 + Sint • (4.25) 

Si (S2) is the effective string field theory action of theory 1 (2). 5i n t is controlled by the 
double-trace angular parameter 9 (see eq. (^j|) ) and includes the interactions generated by 
the double-trace deformation. There are interactions that involve string fields of theory 1 
or 2 separately, but also interactions coupling string fields of different theories (the non- 
vanishing correlation function ((O1O2)) in section || is an example of such interactions). 
5i n t is expected to be a non-local action. Aspects of spacetime non-locality in NLST in the 



context of the AdS/CFT correspondence have been discussed in [47, |4q ]. 

The higher dimensional generalization of this action in AdS/CFT will be a multi- 
gravity theory (when the gravity approximation is justified). 1/N effects will generate 
non-trivial interaction terms. One class of such interaction terms includes a potential for 
the gravitons, in particular a mass term for a linear combination of the gravitons. Such 



effective actions remind of the non-linear bi-gravity actions of 36]. 
5. Discussion 

5.1 Summary and similarities with AdS/CFT examples 



In this paper we found, by generalizing the analysis of j23|, 24 1, that under certain condi- 
tions the product of k large N matrix models coupled together by a multi-trace operator 
admits new double scaling limits, which define holographically k c < 1 non-critical string 
theories interacting with each other in a non-trivial manner. These cases provide one- and 
two-dimensional illustrations of the interacting (multi)-string theories proposed in [Q, §]. 
Specific examples of this phenomenon were given in the main text for k = 2 and in appendix 
[B] for k = 3. At least two conditions had to be met in order to obtain well-defined double 
scaling limits with a non-vanishing multi-trace coupling: (a) the single-trace operators that 
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participate in the deformation should share the same scaling properties, and (b) we should 
fine-tune the multi-trace couplings to a special set of values. In the k = 2 examples of 
section || this special set comprised of a one-parameter family of double-trace couplings 
where the modified double scaling limit was giving rise to a modified partition function, 
which is related to the original product of two one-matrix model partition functions by a 
certain Laplace transform (see eq. ( |1 . 1| ) ) . 

We argued that these deformed large-iV matrix models provide the holographic def- 
inition of an interacting product of c < 1 non-critical string theories. The interaction is 
mediated at the level of the worldsheet via non-local interactions, which induce tiny worm- 
holes connecting worldsheets of the same or different string theories in a manner specified 
explicitly by the exact matrix model free energy. This non-local structure, which is evident 
in the correlation functions presented in section ||, is a special example of the non-local 
string theory construction anticipated on general grounds in |jl9| when multi-trace defor- 
mations are present in gauge/string dualities. 

During the discussion of this structure and the possible interpretations of the modified 
matrix model as a modified string theory we had the chance to revisit an old claim of 
Klebanov p3|| , who proposed that a matrix model deformed by a critical double-trace 
deformation defines a string theory with a wrong branch tachyon condensate. We argued 
that one can make sense of this claim at tree- level as a change of the bulk/boundary 
dictionary, or equivalently as a re-interpretation of the standard minimal string free energy 
in terms of the dual cosmological constant. This observation is analogous to the tree- 
level interpretation of double-trace deformations in the AdS / CFT correspondence as mixed 
boundary conditions for the dual fields in AdS. We can check directly in the matrix model 
case that this picture cannot be extended beyond tree-level, where one should think of the 
modified string theories as bona fide NLSTs. 

The matrix models in this paper and their dual string theories are interesting because 
they offer a set of clean, solvable examples of the general situation outlined in the beginning 
of the introduction: the holographic duality between a product of gauge theories deformed 
by a multi-trace deformation that preserves the product gauge group and string theory 
on a union of spaces. In that sense, the matrix model examples are useful precursors 
of analogous higher-dimensional cases in the AdS/CFT correspondence. The latter are 
interesting because they give us a rare chance to discuss the properties of massive, multi- 
graviton theories in a setting with a UV completion. The higher-dimensional AdS/CFT 
setups are the main focus of a companion paper [2^] . 

Treating the matrix models as a playground for the corresponding AdS/CFT cases is 
further justified by the many similarities that they exhibit. We can summarize the basic 
parallels between the double-trace deformed matrix models and their higher dimensional 
CFT analogs with the following three points. 

Double-scaling limits vs fixed points of RG equations 

The double scaling limits in section [2] zoom around points in parameter space where the 
matrix models exhibit a critical behavior. At special values of the double-trace parameters, 
e.g. when equation ( 2.25| ) holds, a new critical behavior arises, but as we move away from 
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this special submanifold in coupling space either we recover the critical behavior of the 
original undeformed theory, or the theory enters into a branched polymer phase where 
no sensible critical behavior exists. As a necessary condition for the existence of the new 
critical behavior, the scaling properties of the single-trace operators participating in the 
deformation have to be the same. 

A similar picture of critical double-trace deformations arises in higher dimensional 
QFTs p2[ . Let us consider a general double-trace deformation of a d-dimensional CFT of 



the form 



5S 



J d d x [ 5 ii(Oi) 2 +ff2 2 (0 2 ) 2 + 2g 12 1 2 ] , (5.1 



where Oi, 2 are single-trace operators in two theories (1 and 2 respectively) with scal- 
ing dimensions Ai, A2. The deformation ( |5.lD is the direct analog of the matrix model 
deformation ([^) . By analyzing the one- loop /3-functions of the double-trace couplings in 
conformal perturbation theory at leading order in 1/N one finds new fixed points away 
from the origin when the scaling dimensions are either equal or are related by the equation 
Ai = d— A2. We encountered a similar condition on the scaling properties of the deforming 
operators of the matrix models in the previous paragraph. When new fixed points exist, 
there is a one-parameter family of them, which is the analog of the one-parameter family of 



critical behaviors captured by the modified free energy (1.1). Away from the fixed points 
the RG flow drives us back towards the undeformed theory at the origin, or away from 
the origin towards a region where conformal perturbation theory breaks down. In the first 
case, the double-trace deformation acts as an irrelevant operator around the origin - this 
is the analog of recovering the standard double scaling limit in the matrix model. The 
second case is more like the branched polymer phase in the matrix model. 

Beyond tree-level, the RG equations in the higher dimensional cases receive 1/N cor- 
rections which displace the critical circle and introduce a coupling between single-trace, 
double-trace and other multi-trace couplings. At the same time, 1/N corrections in the 
bulk produce an effective potential for the gravitons. In particular, they lift the mass of a 
linear combination of the gravitons and make the dual multi-gravity theory massive. The 
structure of the RG equations on the boundary predicts that the space will remain a union 
of AdS spaces only if we fine-tune the double-trace couplings to a special set of values. 
Otherwise, loop effects will backreact to the background. In those cases, where we can 
trust the perturbation theory and the RG flow on the boundary drives the theory back 
towards the undeformed point, we expect the dual gravity description of the IR physics 
to be captured by a decoupled union of AdS spaces. In this way, the theory develops dy- 
namically a region of spacetime where the mass of the gravitons is washed away. A more 



detailed discussion of these issues can be found in [22|. 



The bulk/boundary correspondence at tree-level 

Another similarity between the matrix model and the higher dimensional CFT cases con- 
cerns the tree-level interpretation of the deformation on the dual string theory or gravity 
side. In both cases the partition function of the deformed theory at the new critical points 
is related to that of the original theory by a Legendre transform. In the bulk, we still have 
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a product of standard minimal string theories or a product of supergravity theories on a 
union of AdS spaces, however, the bulk/boundary dictionary is now modified. In the ma- 
trix model case, we re-interpret the standard amplitudes of minimal strings as a function 
of the dual cosmological constant. In the AdS case, we put mixed boundary conditions on 
the dual fields. 



The NLST structure 

Finally, in both cases and to all orders in the 1/N expansion, the theory in the bulk is an 
NLST along the lines of |19| ]. In the matrix model case, we have an expression exact to 
all orders in perturbation theory for the modified free energy, which is given by a Laplace 
transform. This allows for an explicit derivation of the NLST rules in this case. In the 
higher dimensional AdS/CFT cases, we do not have the luxury of such exact expressions 
but a qualitatively similar structure is anticipated. 



5.2 Open problems and possible extensions 

We would like to close with a short list of open problems. An important issue is the 
non-perturbative stability of the modified string theories in this paper and generalizations 
thereof. The minimal (p, q) bosonic string theories have well-known problems at the non- 
perturbative level, hence Laplace transforms of the form (2.27), ( 2.29Q are at best good 



expressions in a perturbative expansion. It would be nice to re-examine the effects of 
multi-trace deformations in examples which are well-defined non-perturbatively. Perhaps, 
the unitary matrix models of |4(| are a good starting point for such an exercise. 

In this paper we discussed exclusively what happens to closed strings when we couple 
two or more non-critical string theories together. It would be interesting to explore the 



effect of this coupling also to open string sectors. For instance, FZZT branes [49, |5(| 
are represented in the matrix model by the microscopic loop operator logdet(x — <3?). It 
would be interesting to analyze how the multi-trace deformations affect amplitudes that 
involve FZZT branes, i.e. amplitudes that involve a number of insertions of the determinant 
operator det(x — $), and what this implies for the worldsheet theory. 

Another interesting problem is to analyze the pattern of interactions in the case of 
a large number of string theories coupled by non-local interactions. A simple example 
along these lines is to consider M c = matrix models coupled pairwise with double-trace 
interactions in a circular fashion. This is reminiscent of the matrix model couplings that 
define the c = 1 matrix theory in terms of c = matrix models, however, the difference 
here is that the nearest-neighbor couplings are of double-trace type. Such constructions 
should involve a limit where N 3> M 3> 1. 

One of the main motivations behind this work is cosmology and the potential for in- 
teresting time-dependent solutions in multi-gravity theories. One could attempt to analyze 



such solutions in two-dimensional toy examples extending the ideas in [51] to systems of 
coupled c = 1 string theories. 
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Appendices 

A. Products of matrix quantum mechanics models 

In this appendix we discuss double scaling limits in a double-trace deformation of a product 
of two matrix quantum mechanics theories. Double scaling limits in double-trace deformed 
matrix quantum mechanics theories were discussed originally in [52, M]. 
For us the partition function of interest is 1 



J D*!(t)£>* 2 (t) e -^/o^^[ Tr (|*?+l*?- A i*!)+ T >-(|*i+l*i- A ^l)] x 

e -/ 2 ^dt[S 11 (Tr<I.3 ) 2 +2si2Tr$ 3 Tr$ 3 +922(Tr$ 3 ) 2 ] ^ 



The matrix model lives on a compact one-dimensional spacetime with radius R. The 
standard double scaling limit of this model at gn = g 22 = gn = describes holographically 
the decoupled product of two c = 1 non-critical string theories. 

As in the zero-dimensional matrix model case of section |2[ it is convenient to re-exress 
the double-trace deformation as 

g u (Tr $?) 2 + g 22 (Tr $f) 2 + 2g 12 Tr Tr $ 3 2 = 

= ri(CTr$? + STr<S> 3 2 ) 2 + r 2 (-STi$l + CTi$ 3 2 ) 2 , (A.2) 
where again we will use the notation C = cos 9 and S = sin 8. The same relations between 



gij and r\,r 2 and 9 hold as in subsection 2.1 



Following pj], it will be useful to perform the following Fourier transforms 

r 2nR f 27TR nt 

Pi = dt Tr $f , Q „ = V2 / dt cos — Tr , 

Jo Jo & 

r 2nR nt 

S i<n =V2 sin — Tr$f, » = 1,2. (A.3) 

Jo 

Then we can recast the double-trace interactions as an infinite sum of squares 

,-2-kR 

/ dt{CTi^\ + STi^l) = 
Jo 

(CPi + SP 2 ) 2 + Y, [(CCi, n + SC 2 , n ) 2 + (C5i, n + S,S 2 , n ) 2 ] I , (A.4a) 



27T.R 



n=l 



For simplicity, we will set the ranks of the two matrix models iVi and N2 to be equal. This will not 
affect the final results in any significant way. 
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p2itR 

/ dt(-STr$l + CTr$l 
Jo 

-SPi + CP 2 f + £ [(-SCi,„ + CC 2 ,„) 2 + (-SSi,„ + CS 2 ,„) 2 ] \ (A.4b) 



>l) 2 



27T.R 



n=l 



and proceed following the strategy of section ^. 

First, we distinguish between two cases: det g = r\r 2 7^ and detg = r\r 2 = 0. The 
former case is a simpler version of the latter, hence this is the one we will discuss mostly. 
Assuming det g ^ 0, we can use repeatedly the identity ( 2.1 Q| ) to recast the partition 
function (A.l) into the form 



/oo 00 
dyodx Y[ dy n dx n dz n dw n e 

-°° n=l 

/ 



VQ+Vn+ Z n 1 X Q+ X n+ W n 



J X 



D$i(t)D$ 2 (t) e -^/o rf W§* 2 +3* 2 -( A i+ c w-sx )^ x 

e Af E^°=i [(Cy n -Sx n )Ci, n +(Syn+Cx n )C2,n+(Cz n -Sw n )Si, n +(Sz n +Cw n )S2, n ] 5) 

In order to proceed further we need additional information about the partition sum of 
a single c = 1 MQM model. The free energy of this theory is 

lo g y D${t)e- N ti* R dt [ Tr (|* 2 +l* 2 -( A +2/o)'f 3 ]-E^i(?/nC„+ 2 „5„) = 
= 2irRN 2 (- ai x + \a 2 x 2 ^j + F {x,N 2 ) + F l {x,y n ,z n ,N 2 ) , 



(A.6a) 



where 



F (x,N 2 ) = RN 2 -a 3 x 2 /\ogx + --- 



Y4( R+ Ti ] ' 



(A.6b) 



F 1 (x,y n ,z n ,N 2 ) = TrRN 2 f2(yn + z 2 n)(bn + b' n (x/logx) n/R + ■■■)+■■■ , (A.6c) 



n=l 



x = c 2 - A - y 



(A.6d) 



In our case, we have to deal with the decoupled product of two such theories, so we 



set 



y = c 2 - (Ai + Cy - Sx ) , x = c 2 - (A 2 + Sy + Cx ) , 
which can be inverted to 

y = Ai - Cy - Sx , x = A 2 + Sy - Cx . 

In these expressions Ai and A 2 are still defined as in Q2.15a ), ( 2.15b[ ), ( |2.16| ). 
Plugging this information back into (A. 5) we obtain 



(A.7) 
(A.8) 



dxdy JJ dy n dx n dz n dw n e 

-°° n=l 
D 2TTRN 2 (-a 1 x+^a2X 2 -a 1 y+^a 2 y 2 )+Fo(x,N 2 )+Fo(y,N 2 ) x 

F 1 (y,Cy n -Sx n ,Cz n -Sw n ,N 2 )+Fi(x,Sy n +Cx n ,Sz n +Cw n ,N 2 ) 



Vn+'n 1 + 1 (^l-Cy-Sx) 2 , (A 2 + Sy-Cx) 2 



(A.9) 
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The next step is to diagonalize the quadratic terms in the auxiliary parameters and de- 
pending on the sign of eigenvalues take a suitable double scaling limit. 

First, let us consider the quadratic part associated with the auxiliary parameters x,y 

S 2 (x,y) = [(2a 2 det g + g u )x 2 + (2a 2 det g + g 22 )y 2 - 2g 12 xy] • (A.10) 

To diagonalize this expression, we rotate to the variables x± defined by the linear trans- 
formation 

y = U+x + + U y x- , x = U+x + + U~x_ (ATI) 
with f7± and U± the same as Uf- and U± in eq. fl2T9|) . Then we get 

S 2 (x,y) = -N 2 {Mlx\ + M 2 xl) (A.12) 

with 

M| = 2vr Rm\ (A. 13) 

and m± as in eq. (2.21b). 

The corresponding masses for the higher auxiliary modes y n , z n ,x n ,w n are 

Z n , Vn ■ m ly = + b n) > (A.14a) 

x n , w n : ml w = -^^~ + b„\ . (A.14b) 
Since b n < a 2 for all n > 1 [24|, we deduce the inequalities 

™*y>--Z- -+ 2 °2 , mL>-— (- + 2a 2 ) , (A.15) 



2 VP 



which, say assuming r\ > r 2 , imply 



m% > M 2 + , m 2 xw > M 2 _ . (A.16) 

Consequently, the auxiliary parameters z n ,y n ,x n ,w n are massive and can be integrated 
out as long as M± > 0. In that case, the partition function (|A.9| ) takes the simpler form 

dx+dx- e nRN2 ( E+x ++ E ~ x ~) e - N2 ( m + x ++ m - x -) x 

e Fo{U+x + +UxX-,N 2 )+F {U+x++UyX-,N 2 ) ^ (A.17a) 

where 

* = (-2, - ^ + ?f) V} - (*, + ^ + 0* . <A.17b, 

If both m± are positive, we can further integrate out both x± to obtain the direct 
product of two c = 1 partition functions 
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If one of the masses squared is zero, = 0, and the other positive, we can integrate out 
x_ to obtain a Legendre transform formula analogous to eq. ( 2.27| ) 



„ Ar2rl+ *b( U+x++U-^—Z£,N 2 )+F { U+x++U-^A,N 2 

dx+ e vRN Ex +e V J \ y + y ^^' J. (A.19) 

Finally, when both m± = 0, g\2 = 0, gn = 522 = T\ = ^2 = — 20^- 111 this case, we are 
dealing with a decoupled product of two c = 1 string theories that have been deformed 
individually by double-trace deformations to their critical points as in [24|. 

B. Triple intersections 

In the main text we mostly restricted our attention to pairs of matrix models or CFTs 
coupled by a double-trace deformation. In this appendix we consider a triple intersection: 
the product of three theories coupled by a triple-trace operator. As a concrete example, 
we will analyze the case of three 2nd multicritical matrix models. The partition function 
of the theory of interest is 2 

Z= [f[ D^e-^Ef^tl^+^^l-iEL-,^^^^*!)^^)^^) . (B.l) 
J i=i 

The coefficient of the triple-trace interaction has been scaled as iV -1 , so that the whole 
interaction term scales as N 2 . By definition, the parameters gtjk are fully symmetric in 
i,3,k = 1,2,3. 

A general multi-trace deformation that involves n single-trace operators can always 
be recast in terms of single-trace interactions by introducing 2n auxiliary field integrations 



1 54]. In the case of the triple-trace deformation (|B.1| ) this amounts to introducing three 



auxiliary pairs of parameters (<Tj, Uj) and writing 

Z = ff[ D^ i da l dv l e -^E?=i[H+^]-^Ef, j ,^i*ifc^fc+El=i«i(«Ti-^*f) 

i=l 
3 „ 

1=1 ^ 
3 

= II / d(Tidxi e Hxl '° l ' N2) ■ (B.2) 



In the first equality the V{ integrations are defined by analytic continuation. In the second 
equality we changed integration variables from Vi to xi 

Vi = N{ Xi -02- Xi) . (B.3) 

The introduction of the auxiliary parameters (<7j, generalizes the trick of eq. (2.10) to 
an arbitrary multi-trace interaction. 

2 Again for simplicity, we take the ranks of all three matrices to be equal. 
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Saddle point expansion 

The saddle point equations for the integral expression (B.2) are 

1 OF ai 

Xi : -ai + a 2 Xi + — h — = 0, 

N 2 oxj N 



Oi : N(x { - c 2 - Xi) - — 



9iiiO~i + 2y^giik°~i°~k + ^ 9ijkO~jO-k 

k^i 



(B.4a) 
. (B.4b) 



We have denoted the saddle point values of (ai,Xi) as (<7j,Xj). 

The expansion of the free energy J-(xi, o~i, N 2 ) around the saddle point values involves 
the second derivatives 

d 2 T 



dxidxj 



(a F 
N * a2 + ~M 



d 2 T 



doidoj 



6 


" ^ , / 


~ ~N 







k^i,k^j 



i+3- 



d 2 T 



dxidaj 



(B.5a) 



(B.5b) 



(B.5c) 



Several simplifications are possible in the large N limit. First, anticipating the double 
scaling limit, where X{ — ► 0, we can drop the ^rj^. terms in ( |B.5a ) as subleading. Then, 
the leading terms in the saddle point equation (B.4a) give 



(B.6) 



and eq. ( |B.5b| ) becomes 

d 2 T 



da id a j 



= — 6ai 


" , / 

Oijdiii 1 







2^ 9ijk 



-6ai/j, 



(B.7) 



With these simplifications the partition function (B.2) becomes after shifting 



Na 2 



Z = f[[ dx % da % ^i'W+'ZLi ^?-|Eti =1 My^ + oM {m 
i=i •* 



where we define the auxiliary field mass matrix as 



M 2 ij = 6 -^+6a 1 f ij . 
a 2 



(B.9) 



The properties of Z depend crucially on whether the matrix A4 2 is positive definite or not. 
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As an illuminating special case, we will consider in detail what happens when the only 
non-vanishing triple-trace parameter is 5123 = g- The eigenvalues of the mass squared 
matrix 



M' 



a 2 

Qaig 



6aig 
\6axg 6015 ^ 



arc 



1 — 6aici2<? 

a.2 



(double eigenvalue) , A H 



1 + Y2a\a2g 

a.2 



1 



1 



(B.10) 
(B.ll) 

ila\a2 ba\a2 

At the lower end of this interval, g = g* , one of the a eigenvectors becomes massless. At 
the upper end, g = g+, two of the a eigenvectors become massless. The critical values g± 
are therefore good candidates for the definition of a new set of double scaling limits. 

Double scaling limits 

We proceed to analyze the exact partition function 



Ai is positive definite inside the interval 



9- <9 <9+ , 9- 



> 9 a 



3 

Z = Y\ / duidxi 



f =1 [Af 2 (-aia: i + ia2x2)+F(x i ,Af 2 )+Ar(a;i+6 9 af-A i )^]- : |<?^i^2<T3 



i=l 
3 



i=l ^ 



TV 2 ) 



(B.13) 



in a double scaling limit where 

N^oo, Ai = c 2 + Xi + Qgaj ^ (B.14) 

with a particular combination, to be specified, kept fixed. 

It will be convenient to begin with the following manipulations: 

(i) Shift the auxiliary parameters cij — > crj + iV(ai — 02X1). 

(ii) Diagonalize the 3x3 mass matrix A4 2 (see eq. (|B.10[) ) with the linear transformation 3 

Xi =J2 U i X 3- ( B - 15 ) 



In the new variables the action T becomes 

3 

i=l 



^(o-i, Xi, N 2 ) = \F(tfXj,N) - NA^ + N 2 a 2 A i U J i X j - Na 2 o- i U ] i X j 
i=i 

\n 2 (i\ [a_(X 2 + Xf) + A + Xf - 6 5 ai o- x o2 - Na 2 {o- 1 U 3 2 X j + a 2 U(X j ) + (cyclic 



2 

' N 



1, j 



ox - Na 2 U{X-\ [a 2 - Na 2 U 3 2 X j ) (03 - Na 2 U^X j ) . 



(B.16) 



3 The coefficients HI are functions of 01,02,51 whose explicit form we will not determine here. 
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With minor re-arranging, we recast ( |B.16| ) as 



3 i 
F(a h X t ,N 2 ) = Y.FiUfX^N) - NA m - -N 2 a 2 2 [\-(X 2 + X 2 ) + A+Xf] + Gfa) 



i=l 



where 



G(a) 



-S j (a)Xj + W k {a)X 3 X k + RV k X i X j X k , 



-NAidi - Qgai((Ji0 2 + 0-20-3 + cri<7 3 ) 



N 



a\a 2 a 3 



(B.17) 



(B.18a) 



S' (a) = X 2 a 2 A i u{ - Na 2 OiU{ + 6g ai a 2 N \{U 3 2 +Ui)a x + (U(+U 3 3 )a 2 + {U{ +U J 2 )a 3 



+6ga 2 



Ula\o 2 + U ] 2 a\a 3 + U(a 2 a 3 



R jk (a) 



RV k = 6gN 2 a 3 2 U\U{U{i . 



(B.18b) 

(B.18c) 
(B.18d) 

We will consider three distinct cases: (a) both A_, A + > 0, (6) A_ = 0, A+ > or (c) 
A_ > 0, A + = 0. When one of the lambdas becomes negative there will be no double-scaling 
limit, so we will not discuss this situation separately. 

The first case, where both lambdas are positive, occurs when g lies inside the interval 
flB~12l ). To cancel the linear term XjS j in ( |B~T7| ) we shift the X variables setting 



where by definition 



Si 



Xi — Zi + Si 



s 2 



N 2 a 2 2 X 



N 2 a 2 2 \. 



5, 



S 3 



N 2 a 2 X, 



Then we scale the parameters appearing in the free energy in the following way 

Z t = N-^Zi , S = N-^Si . 
The latter follows from the scaling 



(B.19) 



(B.20) 



(B.21) 



(B.22) 



(|B.2l| ) allows us to cancel the explicit dependence of the single-matrix model free energies 
F{U\Zj -\-U\ Sj, N 2 ) on N and to re-express them as functions F(U 3 Zj +Ujsj) of zj, Sj only. 
With this scaling the quadratic and cubic terms R jk XjXk, R ljk XiXjXk become subleading 
and can be ignored. We are left with the single-matrix free energies F, quadratic terms in 
Zi and linear and quadratic terms in Sj. After some trivial algebra, we obtain the following 
expression for the partition function 



i=i 



(B.23) 



where the new coefficients Up d l are functions of ai,a 2 ,g and ij are functions of ai,a 2 ,g 
and 5i - the scaling parameters that were defined in ( B.22| ). The last two Gaussian terms 
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in ( |B,23| ) become delta functions in the large N limit, 4 hence we end up with the partition 
function of three decoupled one-matrix models 



3 



Z(t l ,t 2 ,t 3 ) = He F ^ . (B.24) 



8=1 



Now we come to the more interesting case where one of the lambda eigenvalues vanishes. 
Let us consider first the case with A__ = 0, A+ = 1+1 ^i°2<? = JL > q. Eq. ( pTT7[) simplifies 
a bit to become 



F(ai,Xi, N 2 ) = Y^FiUlX^N) - N&m - ^N 2 a 2 X + X 2 + G{<Ji) 

i=l 

+S j {a)Xj + R jk {a)XjX k + R ijk XiXjX k . (B.25) 



This time we shift the X variables by setting 

X x = Z x , X 2 = Z 2 , X 3 = Z 3 + at2 = Z 3 + S 3 (B.26) 

and then we scale in the following way 

Zi = N'^Zi , S 3 = N~ A / 5 s 3 , oi = N 1 ^ . (B.27) 
The function G(a) breaks up in this limit into two pieces scaling in different ways 

G(a) = N 2 ^G(u) - N^AiUi , (B.28) 

where 

G(u) = -6gai(uiu 2 + u 2 u 3 + uiu 3 ) . (B.29) 

The last two terms in ( |B.25| ) are again subleading in the large N limit and can be ignored. 
Hence, we are left at this stage with the action 



F(ui, Zi , N 2 ) = F(Ui Zj +Ufs 3 ) + N 2 / 5 G(u) - N*' s Ai 



iUi 
i=l 

+N-^{ Zl S 1 + z 2 S 2 ) - l -N 2 ^a 2 \ + z 2 + ~N 2 / 5 a 2 2 \ + s 2 3 . (B.30) 



The z 2 term in (B.30) will contribute a delta function and the z 3 integral will localize at 
z 3 = 0. This simplifies the action J- a bit further to 

3 

F{u i ,z 1 ,z 2 ,N 2 ) = J2F(Ulzi+Uiz 2 +Uis 3 )+N 2 / 5 G{u)-N (i / 5 A i u i 

i=i 

+N-^(z 1 S 1 + z 2 S 2 ) + l -N 2 l\ 2 X + s 2 . (B.31) 



4 We define the large N limit of the last Gaussian factor by analytic continuation when d' are negative. 
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The term N 4//5 (ziS' 1 + Z2S 2 ) has a ^-independent piece 

N 6 / 5 a 2 A i (Ulz 1 +Ufz 2 ) (B.32) 

which we will require to stay finite in the large N limit. This can be achieved with the 
scaling 

WfA, = N- & / 5 5 a , a = 1, 2 . (B.33) 
The remaining combination of Aj's will be scaled in the standard way 

UfAi = N~ 4 ^8 . (B.34) 

Hence, after some algebra T reduces to an expression of the form 

3 § 

T{u l ,z 1 ,z 2 ,N 2 ) = Y F{Ulz l +Ufz2+Uf^—+UtM j u j ) + a 2 5 i {U}z l +U 2 z 2 ) 
z — ' a 2 \ ' 



i=l 



+ 



+iV 2 / 5 [u i {)C li z 1 + IC 2l z 2 ) + ChijU, + T^iH , (B.35) 

where the constants J\f 3 /C m , D 3 \ V 1 (for i, j = 1,2,3, a = 1,2) depend only a\,a 2 . 

To obtain the final result, we diagonalize the quadratic term C ll UiUj and use the large 
./V limit to localize the U{ integrals at U{ = 0. This kills all the contributions which are 
linear and homogeneous in U{. Then, defining the new scaling parameters 

t 1 = a 2 5iU} , P = a 2 5iUf , t 3 = —— (B.36) 

a 2 \ + 

and renaming for aesthetic reasons z%, z 2 as t\,t 2 we get the double-scaled partition function 

/oo 
-00 

This expression, which generalizes eq. (2.27), is a double Laplace transform of the original 



partition function. A double integration is natural, since it is a double eigenvalue (i.e. A_) 
that we tune to zero. 

Similar manipulations can be performed in the last case of interest: A_ = ^j- > 0, 
A+ = 0. We will skip the gory details and present the final, double-scaled expression for 
the partition function 

/oo 
-00 

C. String susceptibility exponents in the continuum approach 



The double scaling limit ( p. 26 ) and the resulting partition function ( 2.27| ) give the holo 



graphic definition of two interacting minimal (2,3) bosonic string theories. In this appendix, 



we present evidence for the tree-level interpretation of these theories in subsection 4.2.1 
by studying the scaling properties of the deformed sphere partition function using the 
continuum formalism of the minimal strings. 



- 39 - 



For completeness, and in order to set the notation, let us begin by recalling a few well- 
known facts about the one-matrix model case. For the main example in this appendix, the 
double scaling limit of the partition function of the 2nd multicritical matrix model 

Z = I e -ATTr[i^+A* 4 ] (ai) 

one finds the free energy 

F(t) = -\t^ - 1 log* + ^U- 5 / 2 + 0(t- w ) (C.2) 

where t is the double scaling parameter. 

When we add a double-trace deformation to obtain the partition function 

Z = j D<5> e -NTr[^ + X^}-g(Tr^ (C . 3) 

and tune g to the critical value — ^ , a new double-scaling limit is possible giving the free 
energy 

This expression can be deduced from ( |C.2j ) and the Laplace transform 

/•oo 

F(t) = log / dt e tl+F{t) . (C.5) 



The string susceptibility exponents 7 are defined at any genus g in terms of the expo- 
nents of the scaling parameters at each order in the expansion of the free energy 

F{t) = ... + # i«*/ 2 + ... , X = 2 - 2g . (C.6) 

The above double scaling limits exhibit different string susceptibility exponents: (C.2) gives 



7 



whereas (CM) gives 7 = ^ 



3" 



There is a simple way to determine the critical exponents 7 from the continuum formu- 
lation of the dual string theory. For (p, q) minimal models coupled to gravity the Liouville 
interaction takes the form 

5SLiouviiie = A* d 2 z O min e a+<t> , a + = - P+ l— (C.7) 
J V^PQ 

where O m i n is the matter primary field with the lowest dimension 

h-min = — • (C-8) 

The Liouville path integral in the partition function can be computed by separating the 
zero mode <j)Q and performing the relevant integral 



r # e^' 2 -» eQ+ '° = — r ( Qt) 
J -00 «+ \ 2a +J 



Qx 

H 2a + . (C.9) 
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In this expression Q = y/2%i= is the linear dilaton slope. Identifying the Liouville in- 



teraction constant \x with the matrix model scaling parameter t we deduce the critical 
exponent 

7 = 2 + —. (CIO) 

For example, when (p, q) = (2, 3) (corresponding to the 2nd multicritical matrix model) we 
have Q = ^ and a + = — and therefore 

7 = = -- (C.ll) 

r p+q-l 2 V ; 

reproducing the above matrix model result. 

For the other double scaling limit (C.4), which involves a critical double-trace defor- 
mation, it has been proposed by ref. |^3[ that one should consider a (2,3) minimal string 
with the wrong branch tachyon in ( |C.7|) , i.e. one should replace a+ — » a_ = — ^§=p 
Indeed, for (p,q) = (2,3) this substitution reproduces the matrix model result 7 = 1/3. 
In subsection |4.2.1| we rephrased this proposal as a transformation from the cosmological 
constant fi to the dual cosmological constant \x. 

We now proceed to apply a similar logic to the two-matrix model case ( |2.27 ). The 
sphere contribution to the free energy 

F(t+,t-) = log dt+ e t + t + +Fx(U+t + +U^)+F 2 (U+t + +U^t.) (ai2) 

J —00 

is the leading contribution in the saddle point approximation. The saddle point value of 
£+ is given implicitly by the following equation 

t + = Uf(U+t+ + f/f £_) 3/2 + U^(U^t + + U 2 tJf 12 . (C.13) 

We should solve this equation for £4. in terms of £4. , £_ and then insert the result into the 
saddle point expression for F 

F(t + , £_) = £+£+ - 2 -{U+t + + U^t^f 2 - 2 -{U+t + + U^t^fl 2 + ■■■ . (C.14) 

It seems difficult to obtain a closed expression for generic £_ , but one can easily deduce 
an expression that involves a perturbative expansion in £_ . Indeed, one can show that the 
sphere contribution to the free energy admits an expansion of the form 

00 _ 

5 — 2n 

F{i + ,t.) = Y,fn(e)t + i *- + •••, (c.15) 

n=0 

where f n (9) are functions of 6 that can be determined. In an effort to reproduce this 
expansion from Liouville theory, we now turn to the continuum formalism of the minimal 
string. 

Before adding the double-trace deformation, the total free energy F(ti,t2) is simply 
the sum of the free energies of the two independent constituent matrix models and the 
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sphere contribution has the following form in terms of Liouville zero mode integrals 

/oo poo Q Q 

#ie*- <lM + / d^-^ ~ h «+ +t 2 " - = tT + tT ■ 
-oo J —oo 

(C.16) 

Rotating to a new mixed basis of coupling constants 

tj = ?7+t + + I7f*_ , t 2 = C/ 2 +t + + u 2 t - ( c - 17 ) 



2 -(l/+J + +C/ 2 -i_)e Q +^ 



we may rewrite the same expression as 

J — oo J —oo 

(C.18) 

To obtain the worldsheet version of the matrix model result ( |C.12[ ) , ( |C.15[ ) we transform 
from t_|_e a+ ^ 1,2 to the dual cosmological constant interactions t+e 01 -^ 1 ' 2 leaving the i_ terms 
invariant. This implies the Liouville interaction 

<5S to tai = f d 2 Zl (u+t+e*-* 1 + f/f t_e Q +^) + J d 2 z 2 (u+t+e^ 2 + U^t-e^ 

(C.19) 

The sphere free energy is then 

F(t + ,t-) ~ /°° dt^* 1 -^**-* 1 -^*-**** 1 + [°° d ^ e Q^-urt + e^.u- t _ e ^2 

J —oo J —oo 

(C.20) 

Expanding the exponentials in powers of t_ and then evaluating the Liouville zero mode 
integrals term by term we obtain 

OO Q + nati OO _ „ 
±- 5 — 2n 

i? (*+'*-)| sphere ~E^)*+ tt e=xXwv (°- 21 ) 

n=0 n=0 

reproducing the matrix model expansion ( |C. 15| ) . We have not attempted to make a precise 
matching of the functions f n (8),h n (9). This would involve going beyond the zero mode 



integrals. We regard the expansion (C.21) as evidence in favor of the picture proposed in 



subsection 4.2.1. 
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